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PREFACE 

This final report describes ot forts on Task 2 which sought to 
determine quantitatively t he optimum number and location of spectral 
bands required to do general rock-type discrimination using thermal 
IR measurements only. This study was performed within the Infrared 
and Optics Division of the Environmental Research Institute of 
Michigan under NASA Contract NAS9-153b2 during the period from 
January - July 1978. Mr. Robert K. Stewart of the Earth Resources 
Program Office at the hvudon 13. Johnson Space Center, Houston, Texas 
was Technical Monitor for this work. 

The work performed under this task represents an effort in which 
special emphasis should be given t o a channel selection while taking 
into account sensor noise, atmospheric absorption as well as the 
available technology for satellite imaging systems in the 1982-83 
time frame. This stuuv would have fallen short of the intended goal 
without the significant contribution ot Dr. Robi rt K. Vincent of Coo- 
spectra Corporat ion. I'ni'or tunat civ . the small number of rock samples 
on which later thermal infrared spectral measurements and a quantita- 
tive mineral analysis had been performed limited this study technically. 

This task was conducted under the guidance ot Mr. Richard R. Logault, 
Head of the Intrared Optics Division. Dr. Quentin A. Holmes was the 
Project Manager for this study. The authors wish to acknowledge the 
substantial contr ibut ions made to this work by Robert Horvath and 
Lest or Wit ter ot ERIM. 
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EXECUTIVE SUMMARY 

Thermal infrared scanners with multiple spectral bands in the 
8-14 pm wavelength region will be important for geological remote 
sensing because this spectral region of the electromagnetic spectrum 
contains unique compositional information about silicate rocks and 
minerals. Of special interest are the emittance minima (reststrahlen 
bands), which occur in igneous rocks at different wavelengths depending 
on silicate rock types. 

A detailed study was performed to determine quantitatively the 
number and location of spectral bands required to do general rock-type 
discrimination from spaceborne sensors using only thermal IR measure- 
ments. The objective of this study was to provide empirical evidence 
for the evaluation of a potential imaging system in the 1982-83 time 
frame. To maximize economic benefits for projected users the following 
specific parameters have been addressed: spectral bandwidths, sensor 

noise, atmospheric absorption/attenuation and temperature variation 
within a scene. 

Beginning with 134 thermal infrared rock spectra in the 8-13 pm 
region from samples with both weathered and polished surfaces, a radia- 
tive transfer model was used to obtain corresponding spectral radiances 
at the top of the atmosphere for these samples under three different 
atmospheric conditions and five sets of random temperatures. Results 
of the current study indicate that both sample temperatures and mineral 
composition parameters (% SiO ? , V^, and M^) can be estimated using 
preprocessed measurements from a spaceborne imaging sensor with three 
1 pm wide spectra l bands . Taking sensor noise into account, bands 
located at 8. 1-9.1, 9.5-10.5, and 11.0-12.0 pm are expected to provide 
imagery with a high signal/noise ratio while permitting discrimination 
of materials with little or no spectral structure in the thermal infrared. 
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Moreover, a pair of fixed linear features (TEMI^ and TEM^) of these 
preprocessed multispectral thermal measurements exist which correlate 
well with only sample temperature (r _> 0.77) and rock-type (r J> 0.72) 
respectively for scenes of geological interest. 

No consideration was given to the deleterious effects of vegeta- 
tion cover other than to suggest a non-linear feature, MD, which would 
exclude such areas from analysis. Adding a fourth spectral band from 
12.0-13.0 pm would complete the spectral coverage of the thermal infra- 
red for other applications and should further improve the ability to 
remotely estimate sample temperatures. The small number of thermal 
infrared rock spectra available from samples whose mineralogical con- 
tent is unequivocally known is unfortunate and should be remedied soon. 
However, this investigation shows that there are strong reasons to 
expect that Earth-orbiting multispectral thermal infrared scanners 
could provide important new information for geologic mapping, while 
simultaneously providing more accurate temperature measurements of 
natural terrain features. 
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INTRODUCTION 


Thermal infrared imaging systems have been the subject of research 
for more than a decade because of the potential they hold for natural 
resource exploration. Remote sensing of rock-tvpos using the thermal 
region depend primarily upon broad spectral exitanee minima of minerals 
in the 8-13 pm wavelength region (Lyon, 1964). While Lyon and Patterson 
(Lyon, et al, I9b9) showed that narrow-bandwidth (A\ l). 1 pm) spectros- 

copy can discriminate successfully between rock-tvpes, a low signal t o 
noise ratio would preclude the extension of this high spectral resolution 
technique t o airborne imaging systems. To overcome this limitation, 
Vincent and Thomson ( 1 c > 7 1 and 1972a) used a ratio image of two 2-3 pm 
wide bands from an airborne thermal infrared scanner to obtain multi- 
spectral data for discriminating between rock types in a quartz sand 
quarry in Oklahoma and over volcanic terrain in Southern California. 
Watson (1°73) developed a heat transfer model for ono-d imens ional 
periodic heating ot a uniformly conductive ground with radiative trans- 
fer to the atmosphere. His model is useful in the interpretation of 
thermal infrared images of rock formations. In particular, he exercised 
this model to predict optimum times for data acquisition in determining 
thermal properties ot terrain. Vincent (1 () 73) demonstrated tlt.u a ratio 
ot two thermal infrared channels could be related to variations in 
silicates, while suppressing temperature variations within the scene. 

More recently Kahle il l, 77) derived an improved heat transfer model ot 
the earth's surface which included what Kahle terms "sensible and latent 
heating". However, because the emphasis was on modeling the heat trans- 
fer of the earth's surface, Kahle ignored t he spectral features ot sur- 
face emissivity and atmospheric emi ss ion/absorpt ion effects on the radi- 
ation received by the remote sensor (i.e., it was assumed that bright- 
ness or radiation temperature determined remotely is equal t o the surface 


3 


temperature). Most of her test area was homogeneous in terms of spec- 
tral omittance (eniissivtty) and she used data collected from an aircraft 
flying at a low altitude. 

An investigation of the improvement in rock type classification 
obtained by incorporating a second thermal infrared band on the Landsat-D 
Thematic Mapper (Holmes and Niiesch, 1978) concluded that more than two 
infrared hands were required to distinguish between variations in general 
rock-type for samples with random temperatures. 

Since a multiband thermal infrared imaging system could play a 
large role in the mapping of rock types from space in the 1980s, it is 
important to determine the number and locations of spectral bands neces- 
sary to discriminate between rock-types successfully in the presence of 
a variety of sample temperatures, and representative effects of the 
earth's atmosphere and sensor noise. The major objective of this inves- 
tigation was to theoretically determine whether rock composition could 
be determined through measurements of infrared radiances in the presence 
of sensor, atmospheric and thermal noise. 



APPROACH 


Before the question of the optimal number, best location and band- 
width of spectral bands for geological remote sensing can be addressed, 
two more basic questions arise: How should non-rock samples be screened 

out and on what basis should an evaluation be made? There is no uni- 
versally accepted set of physical parameters which, being measured 
accurately by a thermal infrared scanner, would be unequivocally use- 
ful for general geological mapping. The procedure we adopted for re- 
lating sensor signals to types of igneous silicate rocks involved the 
following six steps: 

• Separated silicate materials from everything else 

• Used Vincent’s parameters to connect Travis' rock chart 
names to the corresponding mineral composition 

• Connected rock samples with spectral emittance (emissivity) 

• Accounted for temperature of sample 

• Accounted for contribution of atmospheric transmittance 
and path radiance assuming no reflected radiation 

• Accounted for sensor noise 

Having related sensor signals to rock chart names, the criterion used 
for selecting the number and locations of spectral bands was the result- 
ing ability to predict rock-types and estimate sample temperature. 
However, samples used for this study were well characterized and homo- 
geneous. No attempt has been made to account for the mixing of rock- 
types or other terrain elements (such as vegetation) within the instan- 
taneous field of view of a satellite scanner. 

3.1 SCREENING OUT MOST NON-SILICA TERRAIN ELEMENTS 

Prior to applying a procedure for discriminating among rocks with 
different silica content it is desirable to discriminate between silicat 
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rocks and non-silica terrain elements. While silicate rocks exhibit 
broad omittance minima in the 9.0-11.0 pm wavelength region, non-rock 
materials such as vegetation and water have almost no variation in 
spectral omittance in the thermal intrarod region (heeman et al, 1171). 
Carbonate rocks such as limestone and dolomite have very narrow omittance 
minima at 11. A gm and 14.0 ;.m , but otherwise exhibit few emissi'.ity 
variations in the 8-14 ;.m spectral region (Hunt and Salisbury, 19 75) 
making them amenable to the procedure described below. However, phos- 
phates (Hunt and Salisbury, 1^75) and sulphates (Leeman et al, 1971) do 
have spectral omittance minima in the S-14 in region an. can hardly be 
screened out by the following procedure. 

We have developed a non-linear Mineral Discriminant, MD, which 
permits screening out most non-silica terrain elements, except phosphates 
and sulphates as mentioned. MI) is based upon the absence of distinct 
reststrahlen bands in the 9-11 ;.m wavelength region. Since spectral 
exitance is related to mineral composition, we approximated the exitance 
curve of each sample by a parabola fitted through three infrared band 
measurements in the 8-12 ..m range. Initially we used bands at 8. 2-9.0 ,.m, 
9. 0-9. 8 un and 10.5-11.3 i.m; later we used the three spectral bands which 
were found to be "optimum for mapping rock type and temperature". Using 
wavelength in ..m as the horizontal axis and spectral exitance in watts/cm“ 
as the vertical axis, we found that the reciprocal of the focal length, 
f, of this parabola formed an effective discriminant against materials 
with little variation in spectral exitance. Defining our mineral dis- 
criminant, MD, by 

MD = l/f 

we found that all ot our silicate samples had values of MD less than 
0.3, while vegetation and water samples correspond to much larger values 
of MD (indicating a relatively flat omittance curve or little evidence 
of reststrahlen bands). For sulphates with omittance minima in the 
8. 1-9.1 ,.m wavelength region and carbonates with minima in the 11-12 m 

b 


i i i i 4 i k I i I i iii wiki l m 1 i w t 



region, MD would become negative. In addition, a discrimination against 
relatively silica-poor soils and sands with fine particle size such as 
clay rich soils seems feasible (Hunt and Vincent, 1968). 

3.2 USING VINCENT'S PARAMETERS TO CONNECT ROCK CHART NAMES TO MINERAL 

COMPOSITION 

3.2.1 PERCENT OF SILICON DIOXIDE (% Si0 2 > 

In the previous geological applications of thermal infrared remote 
sensing, researchers have sought to exploit the fact that some relation- 
ship exists between the position of the silicate emittance minima 
(reststrahlen bands) and % SiO ? content which could be related to sili- 
cate rock types. Vincent and Thomson (1972b) correlated the "center 
of gravity" of the spectral emittance curves with the ratio of the two 
spectral bands used in their study and showed subsequently a slight 
correlation between the ratios and the SiC^ content of the silicate 
rocks. The reasons for the general correlation between the location 
of the reststrahlen bands and the corresponding Si0 9 conLent of rock 
samples are still not fully understood, although considerable work has 
been done on the spectroscopy of igneous silicate rocks. Our computa- 
tion of the correlation between the wavelength of the emittance minima 
and the percent Si0 0 for twenty-five igneous silicate rocks resulted 
in a correlation of only -0.57. Nevertheless we used % Si0 9 as one of 
our physical parameters because Si0 9 is an important constituent of 
igneous rocks and quantitative measurements of SiC> 2 content by weight 
were available for 97 of the samples. 

3.2.2 V ? M1NERAL0GIC PARAMETER 

There exists no widely accepted quantitative parameter which orders 
igneous silicate rocks according to traditional petrology. Beginning 
with the Travis rock chart (Travis, 1955) Vincent suggested an approach 
using linear combinations of seven minerals (Vincent, 1973). A one- 
dimensional parameter, V^, was created to provide a continuously varying 
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TABLE l 


Terim 


LIMITS OF V y FOR ROCK GROUPS BASED ON 
TRAVIS' IGNEOUS ROCK CHART 


Group 

Equ i granu 1 ar Members 

V ? Range 

A 

Gran it e 

100-84 

B 

Syenite 

84-81 

C 

Nephe 1 ine Sven 1 1 e 

81-80 

I) 

Quartz Monzonite 

80-67 

E 

Monzon i t e 

67-66 

F 

Nephel ine Monzonite 

66-65 

G 

Granod lorite 

65-47 

H 

Quart zdiorlte and Diorite 

47-37 

1 

Gabbro, Diabase and Theralite 

37-18 

J 

Peri dot it e 

1 8-0 
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is defined by the following equation: 


lb 


16 

v B i'"i 

i-i 


whore B = a constant 
o 

Bj * coefficient for the volume percentage ol the 
Ith mineral 


m 


l 


volume percentage ol the itli normative mineral 
in this sample 


(•’) 


The sixteen minerals considered are shown in Tables 2a and 2b. 
Although Travis' e lass i t t cat ion places greatest importance on quart/, 
potass ic feldspar, foldspathoids, pyroxene and olivine with secoinlarv 
importance on the relative amounts of smlic and calcic plagioclaso, a 
stepwise regression indicated that the most significant eight minerals 
for predicting Vincent's infrared ratio were forstorito, anorthito, 
wo l last oni t e, water, magnetite, hematite, albite and quart/ In that 
o rdor . 


1. I SmviKAl DATA SOUROKS 

To obtain the best available data set we began bv manually digi 
tizing l 14 rock spectra in tin* wavelength region t rom S- 1 l pm. These 
spectral omittance measurements include 2b Igneous silicate rocks with 
broken and weathered surfaces (Vincent, 197 0, 7b Igneous rocks tllunt 
and Salisbury, 1974), and 12 metamorphic rocks fllunt and Salisbury, 
197b) . All ol tin* spectra we used t rom Hunt were taken over large 
(representative) surface areas ot samples which had been polished (see 
Append i x) . 
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TABLE 2. MINERALOGIGAL INDEX M 11L 

lb 


!a: m . 

represent volume 

percentages 

of the minerals 

below 

Symbol 

MLnera 1 

Symbol 

Mineral 


m l 

Olivine 

m 9 

Apatite 


m 2 

Wollaston ite 

m io 

Orthopyroxene 

( Magnesium Rich) 

m 3 

Water 

,U 11 

Nepheline 


m 4 

Hemati te 

m 12 

Corundum 


m 5 

Alb ite 

m 13 

Cl inopy roxene 

(Iron Rich) 

m 6 

Quartz 

m l4 

Clinopyroxene 

(Calcium Rich) 

m? 

Olivine 

m 15 

Cl inopy roxene 

(Magnesium Rich) 

m 8 

Ilmeni te 

m 16 

Orthopyroxene 

(Iron Rich) 


2b : B and 

o 

B. have the 

following value 

s (Vincent, 1973) 

B 

o 

= 

1 . 360600 

B 9 

= 

0. 120840 

B i 

= 

-0.002417 

B 10 

= 

0.005732 

b 2 

= 

0.350000 

B U 

= 

-0.006352 

B 3 

= 

-0.013110 

B 12 

= 

0.008251 

B 4 

= 

-0.112930 

B 1 3 

= 

0.189930 

B 5 

- 

-0.002834 

B I9 

= 

-0. 129250 

B „ 

= 

-0.001032 

B 15 

= 

0. 128740 

B 7 

- 

0.022932 

B 16 

= 

-0.011553 

K 

= 

-0.090608 
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3. A SAMPI.K TKMl’KUATURE 

Sample temperature has a greater effect than speetral exitanee on 
thermal Infrared radianee detected hv most mu 1 1 ispeetral scanners. ll 
thermal infrared scanner data are collected on 1 v during periods when 
temperature variations within the scene ot inteiist are small, estima- 
tion ot sample temperature may not he necessary for mapping exitanee 
variations across the scene. In the more likelv case where there are 
significant temperature variations in the scene (several occurrences 
of the same rock type with different surface temperatures) discrimina- 
tion between rock tvpos will be more reliable it some form ot correct ion 
is made on d 'tooted radiance for sample temperatures. As long as a 
given material doesn't undergo a change ot state, the spectral omittance 
(emissivity) is nearly independent of temperature for most common mate- 
rials ami for temperature of the terrestrial environment. Therefore 
the spectral exitanee at any given temperature and wavelength can be 
calculated if the emissivity,' t , is known or vice versa (gravbody 
assumption). However, silicate rocks as selective radiators exhibit 
broad spectral omittance minima in the l )-ll m wavelength region. This 
variety ot spectral structures, which is caused by inter-atomic vibra- 
t ions (Vincent, l l )7S), makes it necessary to have more than one thermal 
infrared band t o estimate, or compensate tor sample temperatures. We 
did not assume that all the samples were at the same temperature. 

Instead we drew a temperature for each sample from a Gaussian tempera- 
ture distribution whose mean value was U)0°K and which had a standard 
deviation ot V’K. Moreover, to enable temperature eltects to be iso- 
lated from omittance minima the drawing ot 1 IA random temperatures was 
done l ive t linos to obtain a data base ot b70 samp It's. 

Emissivity, * ('), is defined as the rat io o! the radiant exitanee 
til a material t o the radiant exitanee ot a blackbodv at the same tempera- 
t ure . 



I.S ATMOS 1MIKR I C TRANSMITTANCE AND PATH RADIANCE 


A thermal infrared radiative transtor model (Amlinp. ot at, 1970) 
was used to transtorm each ot those individual rock spool ra i n t *.■* t ho 
correspond i up spoolral radiance at t ho top ot t ho atmosphoro. Dutinv; 
tlu'so ooir.put at ions wo assumed that in t ho thormal region solar raiii.it ion 
ret looted hy the satnpio is negligible. Figures 1 ami 1? lllustrato the 
spool ra 1 dependence ot atmosphorio t ransnit t t am os ami path radiances 
predicted hy this radial ivo translor modol tor t ho throo atinosphoros 
usod in this stmly. Although o:rono plays an important role in atmos- 
phorio absorption, partionlarlv noar t ho pin wavelength region, 

wator oontont o t tin* atmosplu'fo plays a tnoro dominant rolo in raiiiativo 
t rans t or through tin* atmosphoro in t ho thormal intrarod region. I'lius, 
tomporato wintor (atmosphoro b) roprosonts t ho host viewing condition, 
tho tropioa! atmosphoro (atmosphoro c) t ho poorest , and tin* arot io summor 
atmosphoro (atmosphoro a) roprosonts an intermediate condition tor 
thormal observations tin* this study. 

t.b SFNSOR NOISE 

Sensor signals woro obtained by integral ing spootral railianoos at 
tlu' top ot tlu' atmosphoro ovor oaoli spootral baml. Sonsor no iso was 
simulated by quantizing tho sonsor signals into intervals ot * l J • no iso 
equivalent plioti'iis radianoo dittoroiiv'i' (Holmes, ot al, l l) 78). Wo used 
1.2 v l 0 1 ' photons see” 1 cm - sr" 1 as tho noise equivalent photon radi- 
anoo difference tor oaoli thormal band. 


I 1 


TRANSMITTANCE 












DETERMINATION OF THE NUMBER OF SPECTRAL BANDS 


4.1 AN UPPER BOUND ON PERFORMANCE 

For a given atmosphere and sample temperature, numerical integratio 
of the spectral radiance of each sample as viewed from the top of the 
atmosphere over fifty 0.1 um-wide spectral interval steps in the 8-13 pm 
wavelength region produced a fifty-dimensional sample vector. Such a 
measurement could, in theory, be made by a noise-free "fifty-channel 
thermal infrared scanner". In the calculation of these vectors each 
sample was given five different random temperatures to enable tempera- 
ture effects to be isolated from compositional structure. To us these 
vectors represented a control, i.e., they served to indicate the "best 
achievable" performance when these samples are viewed through the atmos- 
phere without regard to sensor noise or practical limitations on the 
number of spectral bands. A measure of this performance was obtained 
by linear regression of these fifty dimensional vectors for the best 
viewing condition to predict sample temperature, % SiO^, and (M^ 
was only available for Vincent's samples). The resulting r 2 values 
were: for temperature, r 2 = 0.943; for % SiO^, r 2 = 0.739; and for 

V^, r 2 = 0.723. (Note: These r 2 values would not have been unity even 

if the atmosphere had not been included.) Thus our intent was to deter- 
mine the number of spectral bands, hopefully not too large, above which 
the point of diminishing returns is reached. This viewpoint was moti- 
vated by an analogous situation in the visible region of the spectrum 
for Landsat (Wheeler et al, 1976). 

4.2 INDEPENDENCE OF SAMPLE TEMPERATURE AND COMPOSITION 

A principal component analysis on three sets (one for each of the 
atmospheres used in this study) of 670 sample vectors indicated that in 
each case the first three principal components accounted for almost all 
(99.90%) of the data variability. Moreover, in each case a linear 
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combination (hereafter called a linear feature) of the measurements 

could be found which was independent of sample temperature and which 

would correlate quite well with rock-tvpe parameters (% SiO , V_, M ). 

1 7 lo 

Unfortunately, the coefficients of these linear features were different 
for each of the three atmospheric conditions used in this study. 

4.3 PREPROCESSING 

To overcome this limitation and to permit the use of broader 
spectral bands a preprocessing technique which minimized the influence 
of the atmosphere on the ensemble of observations was developed. The 
specific form of preprocessing employed is based upon the role sample 
emissivity and temperature play in radiative transfer through the atmos- 
phere. Explicitly, 

Pn t- h 

Lj j = L(T.) + Lj +0 (reflected radiation) (3) 


where the index i denotes the ith spectral band and the index j denotes 
the j^th picture element in a scene. is the radiance measured by 

the ith band of the sensor when viewing an object with emissivity e 


ij 

and temperature T^ , through an atmosphere of transmittance and path 
radiance (we omit any contributions from reflected radiation). 

Neglecting terms involving the combined effects of emissivity and temp- 
erature, subtraction of the scenic mean from each observation yields a 
value which is proportional to atmospheric transmission times the sum 
of two terms. 


i.e. , L i = L(T) 


Path 


(4) 


so, L i . - L i = (A t + B. AT.)^ 


(5) 


The first term on the right-hand side of Equation 3 is proportional to 
the deviation of the sample emissivity from the mean emissivity of the 


i n 
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scene, while the second term is proportional to the deviation of the 
sample temperature from the mean temperature of the scene. 

Assuming that the standard deviation of the measurements in a 
given spectral band is proportional to atmospheric transmittance, the 
preprocessed measurements are 


Y . , 
ij 



<o i - V 


>/■ 


A, Ac + B . AT . 

ij i j 


( 6 ) 


where 



is defined as the preprocessed signal corresponding to 
the signal in the ith thermal band from sample j. 


I.j. is tlie radiance in the ith thermal hand from sample j, 

L, is the scenic mean for the ith thermal band, 

i 


o 


i 


is the variance of L. in the ith thermal band. 


O 

O' 

n 



is the variance of the signal contributed by the sensor, 

is the deviation of the emissivity in the ith band of the 
jth sample from the mean emissivity of the scene in the 
ith band, 


AT. is the deviation of the temperature of the jth sample 
from the mean temperature of the scene. 


while A. and B, a re constants, 
i i 

Since spectral exitanee is related to composition, we began with 
the idea of approximating the omittance curve ot each sample by an arbi- 
trary choice of six O.S . m wide spectral bands. However, we observed 
that as far as temperature and % SiO, were concerned, three spectral 
bands accounted for the data variability percentages almost as well as 
four, five, or six spectral hands. Thus we focused our attention on 
three spectral bands (8. 2-9.0, 9.0-9. 8, 10.5-11.3 tun) and examined the 
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linear features of this three-dimensional data set for 25 of Vincent's 
samples. The first principal component was found to account for 7 9 • , 
of the variation in this data set. The second principal component was 
virtually independent ol sample temperature (correlation = 0.004) yet 
accounted tor 20. of the variation in the data set. Moreover, when the 
preprocessed signals (Y . ) from Hunt's samples were used in the regression 
developed on principal component one of Vincent's samples we were able 
to predict temperatures for Hunt's samples with a correlation ot 0.71 
and a bias of 2°K. The corresponding predictions of SiO, and values 

using the Y.s from Hunt's samples in conjunction with the second princi- 
pal component from Vincent's samples gave correlations ot -0.77 and -0.59 
respectively. Encouraged by these preliminary results we set about 
determining an optimum set for the three spectral bands. 



‘ 1 
i 


» 

t V 


* 


* 


20 


i 


Jerjm 


5 

DETERMINATION OF THE OPTIMUM THREE SPECTRAL BANDS 

To reduce the computational burden the final locations of the 
three thermal infrared spectral bands were arrived at in stages. 

Because the intent was to determine spectral bands suitable for 
use from space, the atmospheric window from 8-13 um was divided into 
nine 1 pm wide bands as shown in Table 3. It was felt that spectral 
bandwidths of this order would provide imagery with good signal/noise 
ratios as well as provide adequate spectral resolution for discriminating 
mineral composition. Among the eightv-four unique choices of nine spec- 
tral bands taken three at a time, only twenty-seven sets of three spec- 
tral bands were examined. This reduction stemmed from our observation 
that non-overlapping spectral bands were preferable and that it was 
unreasonable to locate the center of the shortest wavelength band above 
9.5 pm or the center of the longest wavelength band below 10.5 pm. (The 
center of the reststrahlen bands usually lies in the spectral region 
from 9. A- 10. 8 pm. ) 

As our standard of comparison we adopted the performance of linear 
regressions developed on preprocessing signals t rom Vincent's samples 
(training) in predicting sample characteristics from the preprocessed 
signals (Y ' ) computed for Hunt's samples (test). Different regressions 
were used for sample temperature and mineral composition. 

Preprocessed signals (Y.) for each of the nine candidate thermal 

bands were prepared using the best viewing conditions (atmosphere b). 

Sensor noise was not added at this stage. Then regressions were used 

to derive linear features of the Y to estimate sample temperature, 

SiO„, V-,, and M,, with each of the twentv-seven sets of three spectral 
2 7 lb 

bands. Each three-band set was judged by how well the linear features 
developed on Vincent's samples (training) performed when applied to 
Hunt's samples (test). The three highest ranked sets were selected 
for more detailed analysis: 
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TABLE 3. CANDIDATE SPECTRAL BANDS CONSIDERED IN THE 
8.0-13.0 ym WAVELENGTH REGION 


Band 

1 

8.1- 9.1 

ym 

Band 

2 

8.5- 9.5 

ym 

Band 

3 

9.0-10.0 

ym 

Band 

4 

9.5-10.5 

ym 

Band 

5 

10.0-11.0 

ym 

Band 

6 

10.5-11.5 

ym 

Band 

7 

11.0-12.0 

ym 

Band 

8 

11.5-12.5 

ym 

Band 

9 

12.0-13.0 

ym 


22 


Lxi inn U- j_lu 


Ll±3L 







Best prediction of % SiO, 7 : 

Band set a: (band 1, band 4, band 7) 

Best prediction of V^: 

Band set B: (band 1, band 3, band 6) 

Best prediction of sample temperature: 

Band set y: (band 3, band 7, band 9) 

The final stage in the determination of the best location of the 

three thermal infrared bands was to test the performance of spectral 
band sets a, B, and y in the presence of sensor noise and variations 
in atmospheric viewing conditions. Sensor noise was simulated by 
quantizing the signal values in each spectral band. Quantized signal 
values for spectral configurations a, B, and y were calculated for all 
13A rock spectra using random sample temperatures to simulate variations 
in ground conditions. The drawing of 13A random temperatures was done 
five times to obtain a quantized data base of 670 samples. Three such 

quantized data bases were prepared — one using atmosphere a, one using 

atmosphere b, and one using atmosphere c. After these quantized signals 
were preprocessed, linear discriminants determined for the Y^s from 
Vincent's samples under the best viewing conditions (atmosphere b) were 
used to estimate sample temperatures and mineral composition from the 
Y^s from first Vincent's and then Hunt's samples under all three viewing 
conditions. 
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RESULTS 

Table 4 shows the regression results (predicted vs. true) for 
Vincent's samples under atmosphere b. Since the training and the test 
data are identical this indicates the best that one can do with three 
1 pm wide thermal bands. For comparison. Table 5 shows the correlations 
which resulted from applying the linear discriminants developed on 
Vincent's samples under atmosphere b (training) to the Y^s calculated 
from Vincent's samples under a different atmospheric condition (test). 
While band set a and 3 retain their position in providing the best 
ability to predict % Si0 9 and even in the presence of sensor noise 
and a change in atmospheric conditions, band set y no longer provides 
the best ability to predict sample temperatures. Table 6 gives the 
correlations resulting when these same linear discriminants are used 
in conjunction with the Y^s calculated from Hunt's samples under all 
three atmospheric conditions. Here was excluded because it could 

not be determined for Hunt's samples. V., values for Hunt's samples 
were estimated by matching the sample names to categories within 
Travis' rock chart and then using a representative value of V^ (e.g., 
Hunt's quartzdiorite was assigned a V^ value of 45, whereas Hunt's 
diorite was assigned a V^ value of 39) . We note that under atmos- 
phere c, the predictions of % SiO^ and V^ for Hunt's samples (last 
three rows of Table 6) are better than for Vincent's samples (Table 5). 
This is probably caused by the polished surfaces of Hunt's samples 
which result in more pronounced reststrahlen bands than the weathered 
surfaces of Vincent's samples. 

Considering the results for a theoretical, noise-free "fifty-channel 
thermal infrared scanner" on our data set, the spectral band set a 
(8. 1-9.1, 9.5-10.5, 11.0-12.0 pm) is seen to provide thermal infrared 
measurements from space with good performance for both sample temperature 


25 






thG 1 






2pm 




t 


TABLE 

4. CORRELATION 

BETWEEN TRUE 

AND ESTIMATE! 

VALUE 

OF TEMPERATURE, 


% 

Si0 2 , V ? AND M 16 

FOR VINCENT 

S 25 SAMPLES 

UNDER 

ATMOSPHERE b 



Band Set 

Temp . 

Si0 2 

V _7 



a 

(band 1,4,7) 

0.952 

0.685 

0.744 

0.964 


6 

(band 1,3,6) 

0.940 

0.685 

0.752 

0.981 


Y 

(band 3,7,9) 

0.978 

0.664 

0.750 

0.952 



( 


i* 




Jm 


TABLE 5. CORRECTION BETWEEN TRUE AND ESTIMATED VALUE OK TEMPERATURE, 

% SiO.,, V, AND M,, . (Estimation of Vincent's 2 5 Samples Under 
2 7 16 

Atmosphere c Based on Training Over 25 Samples 
Under Atmosphere b) 


Band Sat Temp. 
a (band 1.4.7) 0.944 
$ (band 1,3,6) 0.931 
Y (band 3,7,9) 0.87S 


Si0 2 

v 7 

M I6 

0.655 

0.696 

0.912 

0.641 

0.715 

0.947 

0.631 

0.698 

0.901 
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TABLE 6. CORRELATION BETWEEN TRUE AND ESTIMATED VALUE OF TEMPERATURE 
Z SiO., AND V-^ (Estimation of Hunt's 109 Samples Based on 

Training Over Vineent's 25 Samples for Atmosphere b) 





Band Set 

Temp . 

SiO,, 

V 7 

Pred iet ing 

Under 

a 

(band 

1,4,7) 

0.774 

0.729 

0.77 1 

Atmosphere 

b 

8 

(band 

1,3.6) 

0.747 

0.711 

0. 593 



Y 

(band 

3,7,9) 

0.861 

0.624 

0.631 

Predict ing 

Under 

a 

(band 

1,4,7) 

0.794 

0.719 

0.7 71 

Atmosphere 

a 

8 

(band 

1,3,6) 

0.767 

0.547 

0. 595 



Y 

(band 

3,7,9) 

0.881 

0.609 

0.627 

Pred ict ing 

Under 

U 

(band 

1,4,7) 

0. 769 

0.670 

0. 740 

Atmosphere 

c 

p 

(band 

1, 1,6) 

0. 747 

0.651 

0.580 




(band 

3,7,9) 

0.8 <6 

0. 5 U) 

0. 585 






and mineral composition under a wide range of atmospheric conditions. 

We developed two fixed linear features of the preprocessed signals for 
spectral band set a. Feature one, which we call "TEMI^", correlates 
excellently with sample temperature while being independent of mineral 
composition. Feature two, which we call "TEMI 9 ", correlates well with 
sample composition and is relatively insensitive to the temperature of 
the sample. The fixed linear features for spectral band set a were: 

TEMI 1 = 0.5225A(Y 1 ) + 0.56253(Y 2 ) + 0.64071(Y 3 ) 

TEMI 2 = 0.74491(Y 1 ) + 0.06440(Y 2 ) - 0.66405(Y 3 ) 

The correlation matrix between these features and the parameters of 
interest for Vincent's 25 samples is shown in Table 7. 

Figures 3 through 6 show scatter plots illustrating results 
achieved with spectral band set a. In each instance linear regres- 
sions computed using atmosphere b were applied unchanged to samples 
observed under atmosphere a. A scatter plot of predicted and actual 

M. , values is shown in Figure 3. Figures 4, 5, and 6 illustrate the 
lo 

correlation between estimated and actual temperature, % Si0 2> and 
respectively for Hunt's samples. The correlation of 0.794 for tempera- 
ture is excellent considering that the relationship used to predict 
these temperatures were developed from observation of Vincent's samples 
under atmosphere b. 



TABLE 7. 


CORRELATION BETWEEN 
TEMPERATURE, / 


TEMI /TEMI., AND THE PARAMETERS 
Si0 o , V ? AND M 1( . 


(For Vincent’s 25 Samples Only) 



Temp . 

SiO., 

V _7 


TEMI 

TEMI 0 

Temp . 

1.0000 






Si0 2 

-0.0473 

1 . 0000 





V 7 

0.0112 

0.9181 

1.0000 




M 16 

0.0283 

-0.7186 

-0.7678 

1 . 0000 



TEMI 1 

0.9366 

-0.0487 

0.0200 

-0.0808 

1 . 0000 


TEMI 

0.0655 

-0.7087 

-0.7565 

0.9735 

-0.0080 

l . 0000 
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FIGURE 3. CORRELATION BETWEEN TRUE AND PREDICTED 
(Vincent's Samples Only) 
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CONCLUSION 


We emphasize that unfortunately the limited number of thermal 
infrared spectra available, as well as the small number of samples for 
which a quantitative mineral analysis had been performed, prevented a 
more detailed study in the field of general geologic mapping. Never- 
theless, based on these simulations and analyses, the following con- 
clusion was drawn: A thermal infrared imaging system with t hree 1 pm 

wide spectral band s located at 8. 1-9.1 pm, 9.5-10.5 pm and 11.0-12.0 pm 
would provide an opportunity to estimate both sample temperatures and 
mineral parameters (% Si0 o , and M^) from space. 
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GROUPING OF IGNEOUS ROCK SPECTRA (Hunt, 1974) 





X Si0 2 

V-7 



? S10 2 

V-7 

1 . 

Apllte f.ranlte 65 


75.0 

86 

39. 

Quart/ Monzonite 149 

61.8 

73 

2. 

Granite 70 


69.0 

85 

40. 

Mon/onite 154 

63.2 

66 

3. 

Graphic Granite 73 


76.0 

87 

41. 

Monzouiie 173 

65.4 

bb 

4. 

Biotite Granite 76 


67.0 

85 

42. 

Utlte 174 

62.9 

bb 

5. 

Hornblende Granite 

150 

75.9 

85 

43. 

Latite 175 

57.3 

66 

6. 

Granite 162 


64.0 

87 

44. 

Blotite Granodiorite 235 

65.8 

62 

7. 

Granite 245 


72.9 

89 

45. 

Cranodiorlte 64 

61.8 

58 

8. 

Granite 244 


64.0 

87 

46. 

Daclte 56 

65.8 

58 

9. 

Rhyolite 101 


77.8 

88 

47. 

Tridvmite Oacite 399 


45 

10. 

Altered Rhyolite 55 

68.1 

90 

48. 

Hornblende Diorite Porphyry 

404 

40 

11. 

Rhyolite 98 


78.6 

90 

49. 

Hornblende Diorite 152 

48.7 

38 

12. 

Tuff 87 


74.6 

90 

50. 

Diorite 403 


40 

13. 

Green Lapllli Tuff 

89 

69.3 


51. 

Hornblende Diorite 69 

54.2 

40 

14. 

Lapilli Tuff 90 


58.7 


52. 

Diorite Porphyry 44 

56.1 

42 

15. 

Tuff 61 


71.1 


53. 

Hornblende Diorite 240 

56.6 

41 

16. 

Black Obsidian 52 


77.0 


54. 

Quartz Diorite 171 

65.5 

48 

17. 

Brown Obsidian 53 


76.0 


55. 

Andesite 239 

54.1 

41 

18. 

Pumice 62 


70.8 


56. 

Hornblende Andesite 236 

65.9 

39 

19. 

Brown Obsidian 77 


72.2 


57. 

Andesite Porphyry 46 

64.6 

44 

20. 

Vitrophyre 93 


73.1 


58. 

Andesite 121 

58.2 

41 

21. 

Perlite 72 


75.0 


59. 

Andesite 130 

6b. 3 

•44* 

22. 

Syenite 39 


52.2 

83 

60. 

Olivine Gabbro 158 

44.2 

24 

23. 

Syenite 172 


58.5 

83 

61. 

Hypersthene Gabbro 75 

48.2 

30 

24. 

Syenite 170 


60.0 

83 

62. 

Bytownite Calbro 38 

50.9 

34 

25. 

Syenite 178 


66.4 

83 

63. 

Syenite Gaobro 159 

51.3 

3b 

26. 

"Nepheline" Syenite 

83 

59.1 

81 

b4. 

Hornblende Gabbro 132 

45.2 

30 

27. 

Trachyte 109 


66.2 

84 

65. 

Diabase 129 

-.8.2 

28 

28. 

Trachyte 42 


55. 3 

82 

66. 

Diabase 131 

52.9 

30 

29. 

Trachyte 99 


66.9 

84 

67. 

Diabase 2 18 

46.9 

27 

30. 

Trachyte Porphyry 

123 

50.8 

83 

68. 

Diabase 242 

46.1 

27 

31. 

Nepheline Syenite Porphyry 

192 56.0 

81 

69. 

Diabase 155 

48.1 

29 

32. 

Nepheline Syenite 

100 

59.3 

81 

70. 

Olivine Basalt ln6 

46.7 

25 

33. 

Nepheline Syenite 

156 

57.6 

81 

71. 

Flood Basalt 58 

*8. b 

28 

34. 

Phonolite 157 


48.0 

81 

72. 

Peridot Ite 410 


5 

35. 

Phonolite 153 


59.8 

81 

73. 

Peridotite 427 


5 

36. 

Monzonite 233 


53.8 

66 

74. 

Peridotite 118 

42.3 

7 

37. 

Monzonite Porphyry 

234 

62.8 

66 

75. 

Peridotite 317 


5 

38. 

Quartz Monzonite 148 

72.9 

78 

76. 

Peridotite 128 

45.8 

4 
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GROUPING OF METAMORPHIC ROCK SPECTRA (Hunt, 1976) 

77. Pink Marble 360B 

78. Marble 457B 

79. Dolomitic Marble 459B 

80. Dolomitic Marble 458B 

81. Serpentine Marble 460B 

82. Red Quartzite 377B 

83. Purple Quartzite 378B 

84. Red Quartzite 379B 

85. Gray Quartzite 382B 

86. Green Quartzite 376B 

87. Calcium Silicate Hornfcls 471B 

8S. Hornfels 470B 

89. Syenite Gneiss 467B 

90. Albite Gneiss 390B 

91. Felsic Gneiss 310B 

92. Chloritic Gneiss 380B 

93. Augen Gneiss 389B 

94. Hornblende Gneiss 465B 

95. Diorite Gneiss 464B 

96. Sillimanite-Garnet Gneiss 466B 

97. Slate 461B 

98. Chiastolic Slate 462B 

99. Gray Slate 307B 

100. Phyllite 473B 

101. Green Schist 392B 

102. Tremolite Schist 469B 

103. Hornblende Schist 241B 

104. Tourmaline Schist 468B 

105. Hornblende Schist 393B 

106. Chlorite Schist 395B 

107. Crumpled Mica Schist 394B 

108. Graphite Schist 397 
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